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ABSTRACT 

In this study, we compared two different methodologies (standard morphometries and geometric morphometries) 
in order to determine their ability to discriminate the honey bee populations distributed in Iran. In standard morphometries, 
four distances and eleven angles were measured from the forewings of different honey bee populations. A total of 18 
landmarks were utilized for the geometric morphometries analysis. All colonies, described by mean of 40 workers, 
were correctly classified using both methods. In the case of individual wings discrimination, geometric morphometries was 
83.83% successful and standard morphometries was 52.77% successful while the geometric morphometries method is 
much simpler and easier compared to the standard morphometries. Therefore, the discrimination ability of geometric 
morphometries on the honey bee populations was better than standard morphometries. 
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INTRODUCTION 

Two hypotheses have been stated about evolution of honeybees. According to the one hypothesis, honeybees 
evolved in near Ethiopia in Africa and extend over Africa and then Middle East and Europe (Wilson 1971). However the 
second hypothesis explains that, although representatives of most of bee species were original to all the continents, bees 
belonging to the genus Apis were originally to be found only in the Old World, namely southern part of Caspian Sea and 
they were distributed to Africa and Europe (Ruttner 1988). According to Ruttner (1988), there are at least 24 Apis mellifera 
subspecies grouped in three or four evolutionary branches, based on morphometric data. Ruttner el al. (1985) discriminated 
six local populations for A. m. meda, as follows: West and Central Iran (Azerbaijan-Iranian highlands), the subtropical 
coast of the Caspian Sea (Mazandaran), Northeast Iran (Mashad), Southeast Iran (Kerman), Iraq, and Southeast Anatolia, 
from Van Lake to Hatay. 

Morphometries has been used widely for identifying honey bee races, because of high practicability and low 
costs. The first effective procedure, developed in the 1970s, consisted of manually measuring 25 characters of the wings, 
sternites and legs with an ocular micrometer; these measures were then analyzed with multivariate statistics to make the 
identifications (Daly and Baling, 1978). 

The field of geometric morphometries is relatively new (see Bookstein, 1991, Rohlf and Marcus, 1993) and has 
shown very rapid progress over the last few years. As might be expected during a period of rapid development, there can 
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be technical problems in some of the pioneering studies as biologists attempt to apply the new tools. It is only now 
becoming clear how the new techniques should be combined in order to carry out comprehensive analyses of real datasets. 
Bookstein (1996b) gave a list of recommendations for such applications and Bookstein (1996a) gave several 
comprehensive examples of morphometric analyses. However, these accounts do not address some of the types of 
applications that are of particular interest in systematic biology. Bookstein (1994) emphasized that morphometries cannot 
supply homologous shape characters. 

The aim of this study was comparison of effectiveness of geometric morphometries and standard morphometries 
as methods for discriminating honeybee populations based on forewing venation. 

MATERIALS AND METHODS 

A total of 10 colonies were sampled from 10 locations belonging to different geographic regions in Iran 
(Sanandaj: 35 19' N; 46 6' E, Dehgolan: 35 16' N; 47 25' E, Qorve: 35 10' N; 47 48' E, Marivan: 35 31' N; 46 10' E, 
Bane: 35 59' N; 45 53' E, Saqqes: 36 14' N; 46 16' E, Divandare: 35 54' N; 47 01' E, Sarvabad: 35 18' N; 46 22' E, Bijar: 
35 52' N; 47 36' E, Kamyaran: 34 47' N; 46 56') table 1. Honey bees were preserved in preservative solution 
(30 parts distilled water, 15 parts 95% ethanol, 7 parts 38% to 40% formaldehyde and 2 parts acetic acid) to prevent sample 
deformation. The right forewings of 40 worker honeybees per colony were mounted by microscope slides. 
Digital photos were taken from mounted wings using a DP12 camera and SZX12 OLYMPUS stereomicroscope. 

Standard Morphometric Analysis 

The total of 14 morphometric characters were biometrically measured with the tpsdig 2.16 (Rholf, 2010). 

£ 

Morphometric characters included forewing length (FL), fore wing width (FB), cubital index (r) (Figure 1), and eleven 

forewing angles (A4, B4, D7, E9, G18, J10, J16, K19, L13, N23, 026) (Figure 1). The mean values of measurements of 40 
honey bees in each colony were calculated. To discriminate the honey bee populations based on morphometric characters, 
multivariate statistical analyses were performed on measurement values. Statistical analyses were carried out with the 
computer software programs SPSS 17.0 and NtSys Pc. 2. 

Geometric Morphometric Analysis 

Eighteen landmarks on fore wings were digitized on photos by tpsDig 2.16 software (Figure 1). The vein 
junctions were used as homologous points for geometric morphometries (Rohlf, 2010). Later, landmarks were 
superimposed using a generalized procrustes analysis (GPA; Rohlf, 1999). All non-shape variations of these landmarks 
such as orientation (or rotation), scale, and size were removed. Next, digitized landmark outputs were analyzed by 
tpsRelw, tpsReg (Rohlf, 2010), and NtSys Pc. 2. Finally, samples of ten areas were classified with discriminant function 
analysis (DFA) by SPSS ver. 18. A multivariate analysis of variance (MANOVA) was carried out on the landmark data in 
order to compare honey bee populations. 

RESULTS 

Analysis of variance (ANOVA) of morphological characters (standard morphometries) showed that 12 out of 14 
characters (i.e. all except E9 and LI 3) displayed statistically significant differences among honey bee populations 
(P < 0.05). E9 (P= 0.117) and L13 (P=0.08) characters showed no significant differences between populations. Analysis of 
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variance (ANOVA) of morphological characters (geometric morphometries) showed that all variable (lx, ly to 16y) 
displayed statistically significant differences among honey bee populations (P < 0.05). 

All colonies were classified correctly using both the classification functions based on standard morphometries 
(Appendix I). A multivariate analysis of variance (MANOVA) displayed statistically significant differences among honey 
bee populations in standard morphometries (F= 6.74, P=0.00) table 2 and geometric mophometrics (F= 8.03, P=0.00) 
table 3. 

Populations of ten areas based on 32 shape variables of fore wings were classified by Discriminate Analysis 
(DFA) in geometric morohometrics. Cross-validation tests based on DFA correctly classified 83.83% of the colonies 
table 4. DFA showed that honey bee population of Dehgolan was almost differentiated (97.5 %) from other populations 
(Figure 2). Population classification based on standard morphometries (14 variables) showed that DFA correctly classified 
52.77% of the colonies table 5. Statistical analysis of four distances and eleven angles of front wings by DFA 
(standard mophometrics) showed that the honeybee population of Dehgolan was nearly separated (77.5 percent) from other 
areas (figure 2). Therefore, geometric morphometries method could better discriminate honeybee populations than standard 
morphometries. 

Relations of geographical populations were evaluated by UPGMA (unweighted pair group method with arithmetic 
means) method. Results showed that topology of Dendrograms were almost similar. Dendrogram of four distances and 
eleven angles of front wings showed that two populations of Marivan and Sarvabad were similar and two populations of 
Dehgolan and Bijar separated from other populations (Figure 3). Dendrogram of 18 landmarks of Front wing showed that 
two populations of Sanandaj and Qurve were similar and Dehgolan populations completely separated from other 
populations (Figure 4). The both dendrograms based of standard and geometric morphometries, honey bee populations of 
Bane, Qurve and Sanandaj classified each other. 

DISCUSSIONS 

With its wide range of geography, Iran harbors different Apis mellifera meda population adapted to various 
geographical regions extending from east to west and from south to north. Within this huge geography Rutther (1988) 
reported the presence of several Persian honey bee populations. 

Wing shape likely driven by: (1) environmental pressures such as latitude (Alpatov, 1929), altitude 
(Verma et al. 1994; Hepbrun et al. 2000) and climate (Hepburn et al. 2001; Radloff et al. 2005a, b; Tan et al. 2008); (2) 
sexual selection (Radloff et al. 2003a) and (3) abiotic factors such as temperature (Soose 1954) and season 
(Mattu & Verma 1984). Otherwise, wing angles were interrelated with phylogeny rather than phylogeography 
(Diniz-Filho et al., 1999). Farshineh et al. (2007) compared A. mellifera meda populations of Iran (Orumieh, Tbriz, and 
Tehran) with populations in different zones of Turkey (Kiseher and Beypazari) and A. mellifera carnica of northern 
Turkey. They found that honeybee populations in Iran were smaller than honeybee populations in Turkey. 

Morphometries has been used widely for identifying honey bees populations. Morphometric identification 
techniques have improved considerably due to new computational techniques (Francoy et al., 2008). The results of this 
study clearly showed that measurements of size and angle characters can be sufficient to identify or discriminate honey bee 
populations, but the new approach, geometric morphometries, is simple and much more effective than standard 
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morphometries in discrimination or identification of honey bee populations in Iran as well as in the world (Tofilski, 2008; 
Kandemir and Ozkan 2013). 

The results presented here show that geometric morphometries yielded better discrimination of honeybee 
populations (83.83%) than standard morphometries (52.77). Tofilski (2008) demonstrated that geometric morphometries 
(84.9%) yielded better discrimination of three honeybee subspecies (A. m. mellifera, A. m. carnica and A. m. caucasica) 
than standard morphometries (83.8%). Kandemir and Ozkan (2013) compared standard morphometries and geometric 
morphometries ability to discriminate the honey bee populations distributed throughout Turkey. The results showed that 
geometric morphometries was more successful (81.5%) than standard morphometries (70.4%). 

Venation was described either by 11 angles in standard morphometries or by the coordinates of 18 landmarks 
(36 x, y coordinates) in geometric morphometries. According to Bookstein (1991), the use of size-free coordinates instead 
of distances, rotations, or angles leads to more exhaustive descriptions of geometric forms in biology. While standard 
morphometries is restricted to distance and rations of distances, geometric morphometries not only includes these 
measurements indirectly but also allows for wing shape analysis by using the landmark approach. The number of 
characters (36 x, y coordinates) used is thought to be more in geometric morphometries and all the characters were 
size-free; thus, much better discrimination based on wing shape was not surprising. 

CONCLUSIONS 

Geometric morphometries method is much simpler and easier compared to the standard morphometries 
Therefore, the discrimination ability of geometric morphometries on the honey bee populations was better than standard 
morphometries. 

ACKNOWLEDGEMENTS 

The authors thank the Entomology Department, Faculty of Agriculture, University of Kurdistan, Iran, for their 
financial assistance. 

REFERENCES 

1. Alpatov. W. W. (1929), Biometrical studies on variation and races of the honey bee Apis mellifera. L. Quar. Rev. 
Biol. 4; 1-58. 

2. Bookstein. F. L. (1991), Morphometric tools for landmark data: Geometry and Biology. Cambridge University 
Press. 

3. Bookstein, F. L. (1994), Can biometrical shape be a homologous character? Pages 197+ 227 in Homology: 
The hierarchical basis of comparative biology (B. K. Hall, ed.). Academic Press, New York. 

4. Bookstein, F. L. (1996a), Biometrics, biomathematics and the morphometric synthesis. Bull. Math. Biol. 58: 
313+ 365. 

5. Bookstein, F. L. (1996b), Combining the tools of geometric morphometries. Pages 131+ 151 in Advances in 
morphometries (L. F. Marcus, M. Cor ti, A. Loy, G. Naylor, and D. E. Slice, eds.). Plenum, New York. 


Impact Factor (JCC): 1.6864 


Index Copernicus Value (ICY): 3.0 



Comparison of Two Morphometries Methods for Discriminating of Honey Bee {Apis mellifera meda Sk.) Populations in Iran 


65 


6. Daly, H. V. and Balling, S. S. (1978), Identification of Africanized honeybees in the Western Hemisphere by 
discriminant analysis. J. Kans. Entomol. Soc. 51: 857-869. 

7. Diniz-Filho, J. A. , Fuchs, S. and Arias, M. C. (1999), Phylogeographical autocorrelation of phenotypic evolution 
in honey bees ( Apis mellifera F.). Heredity, 83: 671-680. 

8. Farshineh Adi, M. B., Vasfi Gencer, H., Firatli, C. and Bahreini, R. (2007), Morphometric characterization of 
Iranian (Apis mellifera meda). Central Anatolian (Apis mellifera anatoliaca) and Caucasian 
(Apis mellifera caucasica) honeybee populations. J. Apicul. Res. 46(4): 225-231. 

9. Francoy, T. M., Wittmann. D., Drauschke, M., Muller, S., Steinhage, V., Bezerra-Faure, M. A. F., De Jong, D. 
and Goncalves, F. S. (2008), Identification of Africanized honey bees through wing morphometries: two fast and 
efficient procedures. J. Apidol. 39(5): 488-494. 

10. Hepburn, H. R., Radloff ,S. E. and Oghiakhe, S. (2000), Mountain honeybees of Africa. J. Apidol. 31(2): 
205-221. Hepburn, H. R., Radloff, S. E., Verma, S. and Verma, F. R. (2001), Morphometric analysis of Apis 
cerana populations in the southern Himalayan region. J. Apidol. 32(5): 435-447. 

11. Hepburn, H. R., Radloff, S. E., Verma, S. and Verma, F. R. (2001), Morphometric analysis of Apis cerana 
populations in the southern Himalayan region. J. Apidol. 32(5): 435-447. 

12. Kandemir, I. and Ozkan, A. (2013), Comparison of two morphometric methods for discriminating honey bee 
( Apis mellifera F.) populations in Turkey. Turk. J. Zool. 37: 205-210. 

13. Mattu, V. K. and Verma, F. R. (1984), Morphometric studies on the Indian honeybee. Apis cerana indica 
F. - effect of seasonal variations. J. Apidol. 15(1): 63-74. 

14. Radloff, S. E., Hepburn, H. R. and Bangay, F. (2003a), Quantitative analysis of intercolonial and intracolonial 
morphometric variance in honeybees. Apis mellifera and Apis cerana. J. Apidol. 34: 339-352. 

15. Radloff , S. E., Hepburn, H. R., Hepburn, C., Fuchs, S., Otis, G. W., Sein, M. M., Aung, H. F., Pham, H. T., 
Tam, D. Q., Nuru, A. M. and Ken, T. (2005a), Multivariate morphometric analysis of Apis cerana of southern 
mainland Asia. J. Apidol. 36(1): 127-139. 

16. Radloff, S. E., Hepburn, R. and Fuchs, S. (2005b), The morphometric affinities of Apis cerana of the Hindu Kush 
and Himalayan regions of western Asia. J. Apidol. 36(1): 25-30. 

17. Rohlf , F. J. (2010), TpsDig, version 2.16, Software, Department of Ecology and Evolution, State. University New 
York, Stony Brook. Available in http://life.bio.sunvsb.edu/ morph. 

18. Rohlf, F. J. and Marcus, F. F. (1993), A Revolution in Morphometries, Trends Ecol. Evol. 8: 129-132. 

19. Rohlf, F. J. (1999), Shape statistics: Procrustes superimpositions and tangent spaces. J. Classif. 16: 197-223. 

20. Ruttner, F., Pourasghar, D. and Kauhausen, D. (1985), Die Honigbienen des Iran. 2. Apis mellifera meda 
Skorikow, die Persische Biene. J. Apidol. 16: 241-264. 

21. Ruttner, F. (1988), Biogeography and Taxonomy of Honeybees. Springer- Verlag, Berlin Heidelberg, 284 p. 


www.tjprc.org 


editor@tjprc.org 


66 


Javad Nazemi Rafie, Reza Mohamadi & Hossein Teimory 


22. Soose, E. (1954), Effect of temperature on the wing index and chitin colour of the honeybee. Archive fur 
Bienenkunde 31: 49-66. 

23. Tan, K„ Hepburn, H. R„ Radloff, S. E„ Fuchs, S„ Fan, X., Zhang, L. J. and Yang, M. X. (2008), 
Multivariate morphometric analysis of the Apis cerana of China. J. Apidol. 39(3): 343-353. 

24. Tofilski, A. (2008), Using geometric morphometries and standard morphometry to discriminate three honeybee 
subspecies. J. Apidol. 39: 558-563. 

25. Verma, L .R., Mattu, V. K. and Daly, H. V. (1994), Multivariate morphometries of the Indian honeybee in the 
northwest Himalayan region. J. Apidol. 25: 203-223. 

26. WILSON, E. O. (1971), Insect Societies. Harvard University Press, Cambridge. 

APPENDICES 





Figure 1 (a) Forewing Length (FL), Forewing Width (FB), Cubital IndexU); 

(b) Angles of Fore Wing Venation; (c) Eighteen Landmarks 
on Vein Junctions of Fore Wing (A Apis Mellifera 
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Figure 2: Discrimination of Honeybee Populations Based on 
Standard (a) and Geometric Morphometries (b) 



Figure 3: Dendrogram Resulting From a UPGMA Cluster Analysis of Samples From 
Geographical Populations Using Data from the Fore Wing by Standard Morphometries 
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Figure 4: Dendrogram Resulting from a UPGMA Cluster Analysis of Samples from Geographical 
Populations Using Data from the Fore Wing by Geometric Morphometries 


Table 1: List and Code of Collecting Sites of Apis mellifera in Iran 


Population 

Code 

Sample 

Number 

Bane 

Ba 

40 

Bijar 

Bi 

40 

Dehgolan 

De 

40 

Divandareh 

Di 

40 

Qurve 

Q 

40 

Kamyaran 

K 

40 

Marivan 

M 

40 

Saghez 

Sg 

40 

Sarvabad 

Sr 

40 

Sanandaj 

Sn 

40 


Table 2: Multivariate Analyses Based on Standard Morphometries 



Value 

F 

Hypothesis df 

Error df 

Sig 

Wilks lambda 

.137 

6.745 

126 

2863.702 

.000 


Table 3: Multivariate Analyses Based on Geometric Morphometries 



Value 

F 

Hypothesis df 

Error df 

Sig 

Wilks lambda 

.007926 

8.038 

288 

3104.6 

.000 


Table 4: Summary of the Colony Assignments with Respect to Regions 
Based on Standard Morphometries Percent Classifications are 
in Parentheses; N Denotes the Number of Colonies 


Groups 

Ba 

Bi 

De 

Di 

0 

K 

M 

Sg 

Sr 

Sn 

Ba 

12(31.6) 

- 

- 

- 

- 

- 

- 


- 

- 

Bi 

- 

19(47.5) 

- 

- 

- 

- 

- 


- 

- 

De 

- 

- 

31(77.5) 

- 

- 

- 

- 

- 

- 

- 

Di 

- 

- 

- 

28(70) 

- 

- 

- 

- 

- 

- 

Q 

- 

- 

- 

- 

17(43.6) 

- 

- 

- 

- 

- 
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K 

- 

- 

- 

- 

- 

25(62.5) 

- 

- 

- 

- 

M 

- 

- 

- 

- 

- 

- 

15(38.5) 

- 

- 

- 

Sq 

- 

- 

- 

- 

- 

- 

- 

20(50) 

- 

- 

Sr 

- 

- 

- 

- 

- 

- 

- 

- 

21(52.5) 

- 

Sn 

- 

- 

- 

- 

- 

- 

- 

- 

- 

21(52.5) 


Table 5: Summary of the Colony Assignments with Respect to Regions Based on 
Geometric Morphometries Percent Classifications are in Parentheses; 

N Denotes the Number of Colonies 


Groups 

Ba 

Bi 

De 

Di 

Q 

K 

M 

Sg 

Sr 

Sn 

Ba 

35(92.1) 

- 

- 

- 

- 

- 

- 


- 

- 

Bi 

- 

33(82.5) 

- 

- 

- 

- 

- 


- 

- 

De 

- 

- 

39(97.5) 

- 

- 

- 

- 

- 

- 

- 

Di 

- 

- 

- 

36(90) 

- 

- 

- 

- 

- 

- 

Q 

- 

- 

- 

- 

25(64.1) 

- 

- 

- 

- 

- 

K 

- 

- 

- 

- 

- 

33(82.5) 

- 

- 

- 

- 

M 

- 

- 

- 

- 

- 

- 

35(89.7) 

- 

- 

- 

Sg 

- 

- 

- 

- 

- 

- 

- 

34(85) 

- 

- 

Sr 

- 

- 

- 

- 

- 

- 

- 

- 

29(72.5) 

- 

Sn 

- 

- 

- 

- 

- 

- 

- 

- 

- 

33(82.5) 


Appendix I. Classification functions for discrimination of honeybee populations using mean of 40 forewings 
from one colony. The forewings were described by distances and angles between vein junctions. 

Table 6 


Subspecies 

A4 

B4 

D7 

E9 

G18 

J10 

J16 

K19 

L13 

N23 

026 

A. m. meda 

31.015 

106.4 

98.29 

21.115 

93.615 

52.115 

92.663 

78.86705 

14.3021 

89.95199 

34.82199 


Appendix I. (Continued) 


Table 7 


Subspecies 

FWL 

FWW 

CuA 

CuB 

Cl 

A. m. meda 

9.0798 

3.14 

0.5978 

0.25320 

2.3611 
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